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1. Introduction 

From Lattice QCD [jl|] we know that strongly interacting matter of zero net baryon density un- 
dergoes a deconfining and chiral symmetry restoring phase transition from a hadron resonance gas 
at low temperatures to a quark-gluon plasma (QGP) at high temperature. The critical temperature 
of this transition has for several years been quoted [|l|] as T c = 173 ± 15 MeV (with an associated 
critical energy density e c « 0.7 GeV/fm 3 ), but recently this value has been challenged. Since (at 
least for systems with small net baryon density such as those created in heavy-ion collisions at Rel- 
ativistic Heavy Ion Collider (RHIC) energies near midrapidity) the phase transition is not a sharp 
singularity but rather a continuous cross-over where thermodynamic quantities change dramatically 
over a relatively narrow temperature range [Q], the definition of T c is not unique, and one obtains 
slightly different values from different types of observables Depending on the choice of Lattice 
QCD action and extrapolation to the continuum, different groups have recently obtained T c values 
ranging from T c (%yy) = 151 ±3±3 MeV (extracted from the peak in the chiral susceptibility [J3j] ) 
over T c {Xs) ~ T C (P) = 176 ± 3 ±4 MeV (extracted from the strange quark number and Polyakov 
loop susceptibilities [|3|]) to T c = 192 ± 7 ± 4 MeV (extracted from light quark number and Polyakov 
loop susceptibilities [Q]). The differences between the numbers extracted by different groups []3|, |j] 
are larger than the quoted statistical and systematic errors and even exceed the FWHM widths 
AT C ~ 30 — 40 MeV of the various susceptibilities [^] ; they mostly reflect different choices for the 
observables used to set the physical scale and the procedures employed for extrapolating to zero 
lattice spacing [|3],^|]. 

From experimental data collected in yfs = 200A GeV Au+Au collisions at RHIC we know 
that the final hadron abundances from central collisions can be described by a hadron resonance 
gas in a state of approximate chemical equilibrium at r c h em = 163 ±4 MeV, jU# = 24 ±4 MeV, 
and a strangeness saturation factor y s = 0.99 ±0.07 [j|]. The quality of the statistical model fit 
is impressive. Furthermore, the STAR collaboration studied the dependence of the fit parameters 
on the collision centrality and found that neither the temperature r c h em nor the baryon chemical 
potential ji B depend appreciably on the impact parameter [Q] 1 ; only the strangeness suppression 
factor exhibits centrality dependence, beginning at impact parameters > 8 — 9 fm, and drops to 
values around 0.55 in the most peripheral Au+Au collisions [@J. The centrality independence of 
Tchem is in stark contrast to the behavior observed in the same experiment for the kinetic decou- 
pling temperature which is extracted (together with a value for the average radial flow velocity 
(/$) of the fireball at kinetic freeze-out) from the shape of the transverse momentum spectra of 
identified pions, kaons and (anti-)protons T^ n increases significantly with increasing impact 
parameter, from T^ n = 89 ± 12 MeV in the most central to 7ki n = 127 ±13 MeV in the most periph- 
eral collisions. At the same time the average radial flow decreases from (/3) = 0.59 ±0.05 in the 
most central to (/$) = 0.24 ±0.08 in the most peripheral Au+Au collisions, demonstrating a strong 
centrality dependence of the fireball expansion dynamics. 

'The chemical decoupling temperatures extracted from the measured hadron abundance ratios depend somewhat 
on the details of the hadron resonance gas model employed: Ref. [^] gives r c h e m = 163 ±4 MeV, fig = 24 ±4 MeV for 
central Au+Au collisions whereas in Ref. Q one finds r c j, em = 157 ± 6 MeV, fig = 22 ± 4 MeV for the same data. We 
will here use the values from the centrality dependence study presented in Ref. [^] which are consistently on the lower 
end of this range. 
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This characteristic difference in the centrality dependences of the chemical and kinetic decou- 
pling temperatures will be the main focus of this contribution. We will show that the observed 
centrality dependences of the average radial flow velocity and kinetic freeze-out temperature are 
consistent with hydrodynamic behaviour of the fireball medium followed by kinetic decoupling of 
the hadrons from microscopic scattering processes, driven by the collective expansion. We will 
then show that a centrality independent freeze-out temperature is inconsistent with a kinetic de- 
coupling process unless the chemical scattering rates have an extremely (namely almost infinitely) 
strong temperature dependence. We interpret this finding as evidence that chemical decoupling 
of the hadron abundances is driven by a phase transition during which the chemical reaction rates 
decrease precipitously, leaving the system in a chemically frozen-out state at the end of the phase 
transition. Only in this way is it possible to obtain a universal chemical freeze-out temperature that 
is insensitive to the (centrality dependent) collective dynamics and only depends on the thermo- 
dynamic parameters of the phase transition. Obviously, the chemical processes happening during 
the hadronization process itself involve colored degrees of freedom and can thus not be efficiently 
described in hadronic language. We also address the centrality dependence of the strangeness satu- 
ration factor and comment on how our picture also reproduces chemical abundance data measured 
in pp and e + e~ collisions. 

2. Review of proposed phase transition signatures 

In this section we will review several previously published suggestions for QGP phase tran- 
sition signatures in heavy-ion collisions and explain why they either failed or their status has re- 
mained ambiguous. 

2.1 Characteristics of the quark-hadron phase transition 

• Lattice QCD tells us [|I|] that the finite temperature, jUg«0 QCD phase transition is associated 
with a rapid rise of the normalized energy and entropy densities, e/T 4 and s/T 3 , near T c , 
indicating a sudden change of the number of effective degrees of freedom. 

• The square of the speed of sound c 2 s = ^ exhibits a deep minimum near T c , dropping from 
close to the ideal gas value c 2 s = | at T > 2T C to almost 1/10 of that value near T c [jj]] before 
rising again below T c to about half the ideal gas value, as calculated from the hadron reso- 
nance gas model [j7|. (This is followed by a final exponential decrease as T — > which is, 
however, not phenomenologically relevant since the fireball matter freezes out before reach- 
ing such low temperatures). 

• Thermal fluctuations of the net electric charge and baryon number decrease by a factor 2-3 
above T c , as a result of the charge and baryon number being distributed in smaller (fractional) 
units. 

• The transport properties of the matter, such as color conductivity o ab [^] and specific shear 
viscosity r\/s [|I0|], change dramatically at T c . 
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• Finally, chemical equilibration rates, in particular for processes involving the creation and 
destruction of strange quark pairs Qllfl, speed up considerably above T c , due to reduced mass 
thresholds as a result of chiral symmetry restoration. 

All these changes refer to thermal equilibrium properties of the matter created in the collision, 
with the temperature as control parameter. Unfortunately, near T c a few percent rise of the tem- 
perature requires a severalfold increase of the energy density, so these changes set in much more 
slowly when viewed as a function of experimental conditions such as the collision energy, and it 
may be hard to recognize them among other effects. Finally, the medium produced in the collision 
cools through the phase transition quite rapidly, and even if it manages to remain near thermal 
equilibrium in spite of its strong collective dynamics the question remains which of these phase 
transition characteristics survive the expansion and affect the final hadronic state after freeze-out 
in a recognizable fashion. 

2.2 A visit to the graveyard of phase transition signatures 

The first attempt to connect the rapid rise of the entropy density s at T c with experimental 



data was made by van Hove [12] who suggested to plot the average transverse momentum (pj) 
of final state hadrons (as a proxy for the temperature) against their rapidity density dN/dy (as a 
proxy for entropy density JT3|]) - the phase transition should then show up as a "plateau", i.e. as a 



(limited) region where (pr) ~ T remains constant while the entropy density s ~ {\ / nR A ){dN / dy) 



keeps rising with increasing collision energy. However, McLerran et al. soon realized [14] that the 
average transverse momentum (pj) of the final hadrons reflects a combination of random thermal 
and outwardly directed collective motion at freeze-out and is thus not a good proxy for the initial 
fireball temperature. Due to the collective flow component, (pr) depends on the hadron rest mass, 
and since larger initial entropy or energy densities (generated by higher beam energies) increase 
the lifetime of the fireball before freeze-out, this collective flow component keeps increasing as the 
beam energy is raised even when the initial fireball temperature hardly changes inside the phase 
transition. As a result, there is no "plateau" of (pj) vs. dN/dy - the average transverse momentum 
increases monotonically with dN/dy, due to the monotonic increase of the collective radial flow 
with collision energy [jl4j]. This growth is somewhat slowed in the phase transition region where 
the thermal (pj) component remains constant, but that is difficult to observe, and it is even harder 
then to make a convincing case that such a change of slope in {pj) vs. dN/dy is indeed caused by 
a phase transition. 

The speed of sound of the fireball medium represents its "pushing power", i.e. its ability to 
accelerate its collective flow in response to pressure gradients. Near a phase transition it becomes 
small (or even vanishes if the transition is of first order), thus it has been suggested early on (espe- 
cially by W. Greiner, H. Stocker and collaborators [|^]) that phase transitions should be visible in 
the collective flow pattern of the matter created in heavy ion collisions. A specific suggestion made 



by Rischke et al. [ |16| ] within a 1-fluid hydrodynamic model was that the directed flow ("bounce- 
off"), when measured as a function of beam energy, should collapse around Ex^ = 5AGeV due 
to the softness of matter created near the phase coexistence region, but recover again at higher 
beam energies when the fireball is initially in the QGP phase. Subsequent studies within 3-fluid 
hydrodynamics [fL7|], where the two colliding nuclei and the matter created by their collision are 
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described as coupled but separately thermalized fluids, showed that this phase transition signature 
is fragile and exhibits great sensitivity to the non-equilibrium dynamics describing the transfer of 
energy between the three fluids. In a realistic scenario the characteristic minimum in the excitation 
function of the directed flow completely disappears [|l7|]. 

A variation of this theme was proposed in [ |T8| ] whose authors suggested that the softness of the 
Equation of State (EoS) near a phase transition would manifest itself in a non-monotonic elliptic 
flow excitation function, with a minimum at collision energies high enough to initialize the medium 
slightly above the phase transition such that most of the elliptic flow would develop while the matter 
cools through the softest point. The advantage of this suggestion was that elliptic flow would be 
created only after the matter has thermalized, thereby avoiding the very early pre-equilibrium stage 
that proved fatal for the directed flow collapse. However, the elliptic flow minimum turned out 
to be fragile, too, only this time the phase transition signal was killed by non-equilibrium effects 
during the late instead of the early collision stages [19, [ZCl ]; At beam energies below the predicted 
minimum, the elliptic flow would have to be created in the hadron resonance gas phase (whose EoS 
is less stiff than that of the QGP but still much stiff er than the matter near the phase transition), but 
the latter is so viscous that its elliptic flow response to anisotropic pressure gradients is dramatically 
reduced. As a result, the elliptic flow signature never recovers at lower beam energies, thereby 
wiping out the predicted flow minimum around RHIC energies. 

Thus, even though there is now strong evidence from a large variety of measurements that 
the dense matter created in Au+Au collisions at RHIC is initially in the QGP phase [21], no direct 
evidence exists that it indeed passes through a phase transition on its way to the final state of frozen- 
out hadrons. Attempts to capitalize on the predicted change in the net charge and baryon number 
fluctuation spectrum [j|] have not paid off - it seems that the predicted reduction in event-by-event 
fluctuations is largely washed out by the hadronization process. 

Much attention has been paid to non-monotonicities in the excitation functions of several 
heavy-ion observables measured by the NA49 Collaboration in Pb+Pb collisions at low SPS ener- 
gies: a "kink" in the pion production per wounded nucleon, a "step" in the slope parameter of the 
K + transverse momentum spectrum, and a "peak" in the ratio of produced strange to non-strange 
quark-antiquark pairs (for a review see Ref. [22]). The suggested interpretation of these observa- 
tions in terms of a phase transition [^3|, ^2|] is controversial since it relies on a number of unverivied 
simplifying assumptions and is so far not backed up by dynamical model calculations. 

All the above suggestions are based on measuring excitation functions in order to steer the 
medium in a controlled way through the phase transition. The extraction of T c from such studies 
requires a model to quantitatively relate the control parameters y/s or dN/dy to the initial energy 
density (and thus, via the assumption of thermalization), to a temperature). 

We will here propose a different approach that exploits the centrality dependence of observ- 
ables at fixed collision energies. While it has been pointed out repeatedly that the initial energy 
density of the created matter can be varied either by changing the collision energy for fixed col- 
lision system or by changing the size of the colliding nuclei or their impact parameter at fixed 
collision energy, we are not relying primarily on a controlled change of the initial energy density, 
but on the associated variation of the collective expansion rate. In order to explain how we use this, 
let us first discuss the recent controversy relating to the interpretation of chemical freeze-out data 
from relativistic heavy-ion collisions which motivated our work. 
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2.3 The controversy: Chemical reaction kinetics vs. statistical hadronization 

The hadrons emitted in relativistic heavy-ion collisions show thermal characteristics both in 
their abundances and in the shapes of their transverse momentum spectra, but the temperatures 
extracted from yields [||, ||, |4|, 25] ("chemical freeze-out temperature" r c h em ) and from spectra 



VQ, |2g, |27j, [28), |29|] ("kinetic decoupling temperature" T^ n ) differ significantly. The kinetic decou- 
pling temperatures depend on beam energy and, at all beam energies, on system size and colli- 
sion centrality, while for > 10A GeV approximately the same "universal" chemical freeze-out 
temperature is observed in e + e~, pp and A+A collisions at all collision centralities (only the 
strangeness saturation factor y s varies with collision system and centrality) [|(]] . 

This universal value r c h em = 160 — 170 MeV is remarkably close to the critical temperature for 
the quark-hadron transition from Lattice QCD. Furthermore, kinetic simulations of the hadronic 
rescattering stage after QGP hadronization with hadronic cascade codes (RQMD, UrQMD, etc.) 
have shown that hadronic rescattering alters the momentum distributions and resonance populations 
through resonance scattering, thereby cooling the system while keeping it (at least for a while) close 



to local thermal equilibrium [31], while leaving the final stable hadron yields (after resonance 



decays) almost unaffected [ J32| , p3| ]. Hadronic rescattering leads to the loss of a fraction of the 
baryon-antibaryon pairs, but this can be at least partially traced back to the absence of multi-hadron 
collision channels so that detailed balance is violated in baryon-antibaryon annihilation channels 



[34, 35] 



These empirical facts have split the heavy-ion theory community into two camps which offer 



different interpretations of the observations. The philosophy of Camp I is laid out in Refs. [30] 



(second paper) and [36, |33, [?7|] and holds that hadron production is a statistical process associated 
with a phase transition, proceeding through very many different possible microscopic channels 
constrained only by energy, baryon number and strangeness conservation, thereby leading to a 
maximum entropy configuration described by a thermal distribution of hadron yields, with r c h em > 
Hb and y s playing the role of Lagrange multipliers to ensure these conservation law constraints 
while maximizing the entropy. The value of r c h em is not established by inelastic reactions among 
hadrons proceeding until chemical equilibrium is reached - rather, the hadrons are directly "born" 



into a maximum entropy state of apparent chemical equilibrium [|37[|, with the parameter r c h em 
defining the critical energy density at which the hadronization process happens (pO|]b). r c h em 
is thus conceptually different from the kinetic decoupling temperature T\^ n which is the result of 
quasi-elastic rescattering among the hadrons (which also contribute to their collective flow). 



Camp II includes the followers of Refs. Q34[ , p5| , [38| , p9[ ] who hold that chemical freeze-out 
is a kinetic process within the hadronic phase, conceptually equivalent with kinetic freeze-out, the 
only difference being the quantitative values of the corresponding freeze-out temperatures which 
reflect the fact that the inelastic cross sections driving chemical equilibration constitute only a small 
fraction of the total scattering cross section contributing to momentum exchange. The hadrons are 
not born into chemical equilibrium, but driven into such a state kinetically by inelastic multi-hadron 



processes (which, according to Refs. [34, 35, 38, 39], become crucial near T c due to high hadron 



densities) and frozen out by global expansion. Accordingly, r c h em is the "real" temperature at which 
forward and backward chemical reactions last balance each other. (In contrast, for Camp I, there 
are no "backward" reactions involving hadrons in both initial and final states.) 
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Is this more than a philosophical difference of opinions? We think so. Camp II has to cope 
with an intrinsic tension between two observations: The high quality of the thermal model fit to 
the observed hadron yields at RHIC requires sufficient time for inelastic reactions to establish 
chemical equilibrium, whereas the proximity of the fitted chemical freeze-out temperature r c h em to 
the critical temperature T c of the quark-hadron phase transition from Lattice QCD, together with 
the rapid cooling of the fireball by collective expansion, don't provide much of a window for these 
processes to play out. In essence, to make the kinetic chemical equilibration scenario work one 
needs very large scattering rates right near T c which then drop to negligible values just below T c . 
This would be easier to understand if there were a larger gap between T c and r c h em > as suggested 
by the recent upward revision of T c from Lattice QCD advocated in but this problem appears 
serious if the lower T c {Xyy) from Ref. |Q] turns out to be correct. 

In addition, we point out another conceptual problem with the interpretation by Camp II: If 
freeze-out is a kinetic process, it is controlled by the competition between local scattering (mov- 
ing the system towards equilibrium) and global expansion (driving the system out of equilibrium). 
The resulting freeze-out temperature is therefore sensitive to the fireball expansion rate. We show 
that the latter depends on collision centrality. Therefore, any kinetic freeze-out temperature must 
depend on impact parameter. While this is empirically indeed the case for the kinetic decoupling 
temperature T^ m , the chemical freeze-out temperature does not seem to vary with collision central- 
ity. Hence it cannot be the result of a kinetic decoupling process from inelastic hadronic scattering. 



3. Kinetic freeze-out from a hydrodynamically expanding system 

In this section we show that the hydrodynamic model can quantitatively reproduce the ob- 
served centrality dependence of the kinetic decoupling temperature extracted from hadron mo- 
mentum spectra at RHIC. We then show that an analogous centrality dependence of the chemical 
freeze-out temperature cannot be avoided if the hadron yields are similarly controlled by kinetic 
freeze-out from inelastic hadronic rescattering. 



We use our (2+l)-dimensional longitudinally boost-invariant hydrodynamic AZHYDRO [ ]4Q| ] 
with standard initial conditions [f7|] to generate the flow pattern for 200A GeV Au+Au collisions. 
This code has been previously shown to successfully reproduce the measured single particle hadron 
^r-spectra and their elliptic flow (for details see pl|]). Here, however, we modify the freeze-out 
criterium for kinetic decoupling to account for its kinetic nature: Instead of requiring freeze-out on 
a surface of constant energy density edec = 0.075 GeV/fm 3 (corresponding to a fixed temperature 
Tidn = 100 MeV [0]), we define the kinetic freeze-out surface as the set of points where the local 
expansion rate 1 / T exp (x) = d ■ u(x) (u^ (x) being the hydrodynamic flow 4- velocity) becomes equal 
to the local scattering rate 1/T scatt (x) [H|, fj|, JHj]: 



= £, — =£d-u(x). (3.1) 

^scatt ^exp 

Here ^ is a proportionality constant of order unity which we set in a first attempt to £ =0.35, yield- 
ing an average temperature along the freeze-out surface for central Au+Au collisions of (7]dn) — 



115 MeV. Having fixed E, in central collisions, Eq. (3.1) is taken to define the freeze-out surface 



also at other impact parameters. Inside the freeze-out surface the scattering rate exceeds % times 
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the expansion rate, and the matter is thermalized, whereas outside the surface the expansion rate 
exceeds ^ i times the scattering rate so we assume the hadrons there to be decoupled from the 
fluid, streaming freely into the detector. 

The expansion rate d ■ u = 7j_ + V ■ v_|_) + (d t + v^-VjJ /j_ is computed from the hydrody- 
namic output for the transverse flow velocity Vj_(jc) and /j_ = (1 + v 2 ) -1 / 2 . 

The scattering rate in Eq. ( |3~l"| ) involves cross sections and densities of scatterers, hence it is 
particle specific [43|, 44]. Since at RHIC energies hadron production is dominated by pions, we 
assume for simplicity that all hadrons decouple when pions freeze out. It may be somewhat better 
to describe different hadron species by different freeze-out temperatures [44], but each of them 
would exhibit similar impact parameter dependences as the one for pions which we work out here. 

The pion scattering rate we take from the numerical results presented in Ref . [ 44 ] which we 
parametrize as 

1 / T \ 3 45 

-_ = (59.5fm- 1 ) — — . (3.2) 



KieV 

This defines the momentum exchange rate to be used for the calculation of kinetic freeze-out and 
will need to be modified when discussing chemical freeze-out below. 



freeze-out surface for Au+Au 
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Figure 1: Kinetic freeze-out surface Tki n (r) for central (b — 0) 200 A GeV Au+Au collisions, computed from 



Eq. (3.1 ) with % = 0.35 (red solid line) and for a constant freeze-out temperature T^ n = 115 MeV (dotted 
black line). Both surfaces have the same average temperature of (T) = 115 MeV, using the energy density 
as weight function. 



In Fig. H we show the kinetic freeze-out surface for central Au+Au collisions computed from 
0.35 (solid red line) and from the condition Ty^ 



Eq. (3T ) with £ = 0.35 (solid red line) and from the condition 7^, = 115 MeV (dotted black line). 
Both have the same average kinetic freeze-out temperature, but for the kinetic freeze-out criterium 



(3.1) the middle of the fireball freezes out a bit later at lower temperature and larger flow whereas 
the edge decouples earlier at higher temperature and with less flow than the contant-J surface. This 
is caused by the larger expansion rate near the edge of the fireball. 
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Figure 2: Impact parameter dependence of the average kinetic decoupling temperature (T^ n ) (left) and 
average radial flow velocity (vj_) along the freeze-out surface computed from Eq. (3.1) with | = 0.35 and a 
scattering rate ~ T 3 45 as in Eq. (3.2), for 200A GeV Au+Au collisions. 



Figure ^ shows the impact parameter dependence of the average kinetic decoupling temper- 



ature and the associated average radial flow calculated from the kinetic freeze-out criterium (3.1) 
with <^ = 0.35. One sees that central collisions decouple at relatively low temperatures with large 
average radial flow whereas peripheral collisions freeze out earlier when the fireballs are still hotter 
and less radial flow has developed. This is in good qualitative agreement with the STAR data in 
Ref. [||], although there the freeze-out temperatures are generally a bit lower, with slightly larger 
average radial flow velocities than seen in Fig. 0. We adjust for this by fine-tuning the phenomeno- 
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Figure 3: Impact parameter dependence of the average kinetic decoupling temperature (T^ n ) computed 
from hydrodynamics with kinetic freeze-out criterium (3J_) using B, = 0.295, compared with STAR data [|| 
for 200A GeV Au+Au collisions. 



logical parameter £, in Eq. (3.1) to B, = 0.295 (<^ _1 = 3.4). The corresponding freeze-out temper- 
atures are shown as a function of impact parameter b in Figure |3| together with the STAR data. 
Now the agreement is also quantitatively acceptable. We conclude that the measured centrality de- 
pendence of 7ki n can be completely understood in terms of a hydrodynamic model for the fireball 
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expansion coupled to a kinetic freeze-out criterium with realistic temperature dependence of the 
microscopic scattering rate. 

Let us now see whether we can similarly understand chemical freeze-out as a kinetic decou- 
pling process from inelastic hadronic scattering. A few typical processes relevant for chemical 
equilibration are 

K + K< — >K + K, K + Ni — >K + Y,, K + Y< — >K + N, 

Q + x< — >'E + K, K + Ki — >(j) + K, a + K< — >S + 7T, (3.3) 

Q. + N< — >2ti + 3K, N + N< — >5k, N + 3K< — >Q. + 37i. 



The last line shows so-called multi-hadron collision channels which in at least one direction require 
collisions between more than two hadrons. Rates for processes involving n m incoming hadrons are 
proportional to the product of their densities ~ Y^^n^T) where each factor nj(T) grows with 
T at least as T 3 (even much more rapidly for hadrons with masses > T). At low temperatures, 
multi-hadron collision processes as well as collisions between very massive hadrons are therefore 
strongly suppressed. Consequently, particle yields for hadrons requiring collisions of many abun- 
dantly available particles for their production or destruction (such as p,£l,.. .) thus tend to freeze 
out at higher T than particle yields for hadrons whose abundances can be efficiently changed by 
two-body reactions {%,K,§,.. .). 

In an expanding, cooling system, simultaneous freeze-out of all hadron yields at a common 
temperature therefore requires a conspiracy of rates with widely differring T-dependences. Indeed, 
thermal model fits to hadron abundances with a single common temperature are usually not perfect 
[45], and individual fits to subsets of yields measured in lower-energy collisions at the SPS and 
AGS tend to lead to a significant spread of chemical freeze-out temperatures [46]. Only at RHIC 
energies so far the single-temperature chemical equilibrium fit gives an almost perfect fit to the data 

One way to achieve the conspiracy of different chemical equilibration rates that is required for 
a good fit with a single freeze-out temperature is to postulate that at chemical freeze-out all chem- 
ical reactions are completely dominated by multi-hadron collisions and that at any temperature 
below r c hem the medium is so rarefied and so rapidly expanding that even the simplest two-body 
reactions among the most abundantly produced hadrons (such as those listed in the first line of 
Eq. ( |3.3[ )) have essentially stopped. As long as collision channels with widely different tempera- 
ture dependences compete with each other, chemical freeze-out of all hadron species at a single 
temperature appears to impossible. 

Even more importantly, even if it were possible at one fixed impact parameter to arrange for 
common freeze-out of all hadron species in spite of a competition of scattering rates with different 
temperature dependences, such a conspiracy would be impossible to maintain, with the same value 
for the freeze-out temperature r c h em > over the entire impact parameter range. Figure ^ shows the 
centrality dependence of the average chemical freeze-out temperatures along hydrodynamic decou- 
pling surfaces computed with the kinetic freeze-out criterium (3.1), using t, = 0.95 to adjust the 
value of (r c hem) hi central Au+Au collisions to the STAR data [|6|] and exploring different possible 
temperature dependences of the dominant inelastic scattering rate. One sees that approximate im- 
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Figure 4: Impact parameter dependence of the average chemical decoupling temperature (r c t, em ) computed 
from hydrodynamics with kinetic freeze-out criterium (3.1 ) using | = 0.95 and reaction rates with different 
temperature dependences as listed, compared with STAR data [||] for 200A GeV Au+Au collisions. 



pact parameter independence of (r c hem) can only be achieved if all inelastic scattering rates grow 
with T as T" with a power n > 20! 

Basically Figure |] tells us that the observed centrality independence of r c h em requires chem- 
ical freeze-out to happen in a region of parameter space where all chemical reaction rates exhibit 
extremely steep temperature dependence, dropping like a stone as the system cools through the de- 
coupling temperature. It is hard to understand such a behavior within a hadron rescattering picture 
unless one assumes that all relevant chemical reactions involve multi-particle channels involving 
many hadrons. Making such an assumption clearly pushes the hadronic rescattering model towards 
breakdown because its chemical kinetics would essentially be controlled by interactions among 
clusters of particles involving an unspecifiable number of hadrons. It is much more natural to 
associate this kind of behavior with the quark-hadron phase transition where densely spaced and 
strongly interacting quarks and gluons provide the necessary multi-particle clusters, and where the 
dramatic change in number and quality of the effective degrees of freedom within a narrow tem- 
perature interval generates the dramatic temperature dependence of the chemical reaction rates at 
decoupling which seem to be phenomenologically required. 

In such a picture, hadrons are not really well-defined states until after the quark-hadron phase 
transition is complete and, at the same time, chemical reactions among hadrons have ceased. 
Hadrons are thus indeed "born into chemical equilibrium" [37] in a process that can be rightfully 
called "statistical hadronization" [3C, 33]. If hadrons are formed in this fashion, their measured 
abundances provide a window with a direct view of the QCD quark-hadron phase transition. 



4. Conclusions 



We have shown that the observed impact parameter dependence of the average temperature 
and radial flow velocity at kinetic (thermal) freeze-out (i.e. at the point where the hadron mo- 



ll 



Universal chemical freeze-out as a phase transition signature 



Ulrich Heinz 



mentum distributions decouple) can be quantitatively understood as a kinetic decoupling process 
in a hydrodynamically expanding source, with freeze-out being driven by the global expansion 
of the collision fireball. Any such kinetic decoupling process is controlled by the local competi- 
tion between temperature dependent scattering and hydrodynamic expansion rates, and since the 
latter change with impact parameter as a result of the varying initial energy density and size of 
the nuclear collision zone, the resulting average freeze-out temperature is necessarily impact pa- 
rameter dependent. The strength of this impact parameter dependence (i.e. the sensitivity of the 
freeze-out temperature to the fireball expansion rate) is inversely related to the strength of the tem- 
perature dependence of the local scattering rate. To obtain approximate centrality independence of 
the freeze-out temperature, the scattering rate must exhibit an almost infinitely steep temperature 
dependence. 

From this it follows that the observed impact parameter independence of the chemical freeze- 
out temperature in Au+Au collisions at RHIC (i.e. of the temperature where the abundances of sta- 
ble hadron species decouple) cannot be consistently described as the result of a kinetic decoupling 
process from inelastic hadronic interactions. To obtain the necessary extremely steep temperature 
dependence of the inelastic scattering rate (~ T" with n > 20) requires that at the freeze-out point 
all chemical reactions are dominated by multi-hadron interactions involving many more than two 
colliding particles, in which case it seems unlikely that one will ever be able to describe this process 
quantitatively in hadronic language. 

In our opinion the only theoretically consistent interpretation of the STAR data on chemical 
freeze-out is to associate the steepness of the the temperature dependence of chemical equilibration 
rates with a phase transition (in this case the quark-hadron transition) in which the hadrons are pro- 
duced statistically and distributed among different species according to the principle of maximum 
entropy, via a multitude of complicated microscopic channels involving large numbers of strongly 
interacting quarks and gluons. In this sense the hadrons are "born into chemical equilibrium" in an 
environment that is too dilute and expands too rapidly to allow for any further inelastic reactions 
among the hadrons. 

Tidn and r c h em thus stand on conceptually different footings. r c h em is a Lagrange multiplier 
related by the Maximum Entropy Principle to the critical energy density e c for hadronization. Its 
universality in e + e~, pp, and AA collisions of all centralities shows that at e c a phase transition 
occurs. Hadrons are formed during this transition in a statistical process subject to the Principle of 
Maximum Entropy. 

The absence of inelastic hadronic rescattering processes allows the direct measurement of 
T c through r c h em and thus the experimental observation of the phase transition. In this context 
the question arises which of the different definitions of the critical temperature T c from Lattice 
QCD that were mentioned in the Introduction is most closely related to the chemical freeze-out 
temperature r c h em extracted from hadron yield data. It seems unlikely that hadron yields can be 
considered frozen out before the hadrons have more or less recovered their full vacuum masses, and 
this is related to the restoration of the chiral condensate (t/n//) to its vacuum value. We therefore 
suggest that T c {%y V ) ^ should be the LQCD number most closely related to the phenomenological 
value r c hem- This seems to be consistent with the actual values extracted in ^ from LQCD and in 
[g] from hadron yields at RHIC (both are between 150 and 160 MeV). 

The increase of the strangeness saturation factor y s from e + e~ and pp to heavy-ion colli- 
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sions and from peripheral to central Au+Au collisions at RHIC shows that the lifetime of the QGP 
(and thus the time for chemically equilibrating strange with light quarks) is still limited. Only for 
midcentral to central Au+Au collisions y s has sufficient time to saturate. (Qualitatively similar ten- 
dencies are seen in Pb+Pb collisions at lower SPS energies [25].) The primary parton production 
process at the beginning of the collision apparently suppresses the production of strange quarks, 
and it also produces s and s locally in pairs, thereby generating spatial correlations among s and s 
which ensure strangeness conservation locally. In a grand canonical description such correlations 
induce a strangeness suppression factor y s < 1 [BO]. It takes time to diffuse the strange quarks 
over the entire fireball volume to decorrelate them and adjust their abundance to equilibrium val- 
ues. Larger initial energy densities in central Au+Au collision provide more time until the point of 
hadronization at e c ~ 0.7 GeV/fm 3 is reached than peripheral Au+Au or e + e~ and pp collisions. 

We close by pointing out that our conclusions about the nature and origin of r c h em can be put 
to a relatively easy experimental test: It is well known [^, 25] that at low SPS and AGS energies, 
where the net baryon density of the matter created in the collision is much larger than at RHIC, the 
measured chemical decoupling temperatures are well below generally accepted estimates for the 
phase transition temperature, r c h em < T c . In that case the phase transition can not be the origin of 
the observation of chemical equilibrium yields; hadronic chemical reactions must be responsible 
for lowering the chemical freeze-out temperature to values significantly below T c . Since the present 
work has shown that the kinetic decoupling of hadronic chemical reaction rates is influenced by 
the fireball expansion rate, which again depends on collision centrality, we expect to see impact 
parameter dependence o/T c h em whenever its value is measured to be well below T c . This conclusion 
would also apply to RHIC collisions if Lattice QCD would eventually converge to T c values above 
190 MeV as proposed in In this case we would definitely expect r c h em to depend on collision 
centrality and ask for a reanalyzis of chemical decoupling data at RHIC with higher statistics in 
order to unambiguously settle this question. 

It may be possible to re-analyze existing SPS data to confirm or falsify our prediction of 
centrality dependence of r c h em at these energies. If not, this will be a worthwhile point to address 
within the planned low-energy collision program at RHIC. Clarification of this point will be of 
utmost importance for establishing the observed chemical decoupling temperature at RHIC as a 
direct measurement of the critical temperature of the quark-hadron phase transition in QCD. 
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